Staphylococcus aureus and Enterococcus faecium isolated from pigeon droppings (Columba livia) in the external environment close to hospitals

ABSTRACT Background: Domestic pigeons carry pathogens in their droppings, posing a potential public health problem. Methods: The phenotypic and genotypic antimicrobial resistances of Staphylococcus aureus and Enterococcus faecium in the feces of urban pigeons near hospitals with intensive care units were measured. Results: Twenty-nine samples showed Enterococcus growth, whereas one was positive for S. aureus. The S. aureus isolate was sensitive to the antibiotics tested via antibiogram, however resistance genes were identified. E. faecium isolates showed phenotypic resistance to gentamicin, erythromycin, and ciprofloxacin. Conclusions: Antimicrobial profiles harmful to health were demonstrated in bacterial pathogens isolated from the external environment of hospitals.

Domestic pigeons (Columba livia Gmelin, 1789) are the most common animal species in urban environments, where they can find abundant food sources, shelters, and do not have natural predators 1 . Pigeons are often reported as reservoirs of various zoonotic pathogens 2,-5 ; therefore, their presence near hospitals and other patient care settings can be dangerous because of the possibility of pathogen dissemination within the healthcare environment, especially among immunocompromised populations such as intensive care unit patients.
Healthcare-associated infections are frequent adverse events, often resulting in prolonged hospitalization time, treatment failures, post-surgical and post-organ transplantation complications, psychological stress, increased mortality, and increased costs for healthcare systems 6 . Concerning the etiology of nosocomial infections, Enterococcus spp. and Staphylococcus aureus are the most frequently reported gram-positive organisms in epidemiological studies from recent decades. This is mainly due to the increasing prevalence of their antibiotic-resistant variants, such as methicillin-resistant S. aureus and vancomycin-resistant Enterococcus 7 .
This study aimed to isolate and characterize the antimicrobial resistance profiles of S. aureus and E. faecium from pigeon droppings collected within the vicinity of hospitals equipped with intensive care units in Macrorregião Sul do Triângulo Mineiro, Brazil.
For 3 consecutive weeks, pigeon stool samples (n = 96) were obtained from the external areas of 4 hospitals in the south of Triângulo Mineiro (Minas Gerais, Brazil) between July and August 2018. For each facility, 8 weekly collections were made (alternating 4 dry and 4 fresh samples), at an interval of 48 h, avoiding repeated sampling from the same excreta. Samples were inoculated on Baird-Parker and Enterococcosel agar media (MBiolog Diagnosticos, Contagem, MG, Brazil), and cultured at 36.0 °C for 24-48 h. Staphylococcus spp. were identified by Gram staining, colony morphology, and catalase and coagulase tests. Enterococci were identified by Gram staining, 6.5% NaCl broth growth, bile esculin hydrolysis, L-pyrrolidonyl-β-naphthylamide hydrolysis (PYR), and motility tests (MBiolog Diagnosticos). Positive isolates identified from these tests were characterized using a commercial kit for Enterococcus identification (Probac, São Paulo, Brazil) following the manufacturer's instructions. Viable cells were counted using the pour plate method, and the findings were expressed as colony-forming units (CFU/mL).
DNA was extracted and purified using a QIAamp DNA Stool Mini Kit (Qiagen, Venlo, Netherlands). The amount of DNA in the samples was determined using a NanoDrop spectrophotometer.
Quality (integrity) analysis, quantification, library preparation, and sequencing were performed at the LaCTAD laboratory (Campinas, Brazil). Quality was measured using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Absence of degradation peaks or dragging of degraded DNA was observed, and quantification of intact samples was determined using a Qubit 2.0 fluorometer (Thermo Fisher Scientific, Waltham, MA, USA). The genome library was prepared using the Nextera XT DNA Library Preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer's guidelines, and evaluated using Bioanalyzer. Then, these were quantified in Qubit, underwent quantitative PCR with the KAPA Fast Universal kit, and sequentially diluted to prepare the pools. Subsequently, flow-cell clustering was performed, followed by sequencing on HiSeq2500 (Illumina), with an expected coverage of approximately 10 million reads per library. Sequencing reads were processed using 3 tools from the BBTools version 38.79 package (https://sourceforge.net/projects/bbmap/, https://jgi.doe.gov/data-and-tools/bbtools/). Illumina adapters and low-quality bases were removed using BBDuk (k = 23; mink = 9; hdist = 1; ktrim = r tpe tbo; minlen = 135; qtrim = rl; trimq = 15; maq = 25). BBMerge was used to correct errors in overlapping regions between the forward and reverse reads (ecco = t; mix = t; strict = t), while Clumpify was used for another round of sequencing error correction. An assembly was then made with metaSPAdes version 3.14.1, using the parameters: k = 31, 51, 71, 91, and 111. Contigs smaller than 200 bp were removed.
The genome draft underwent contaminant analysis using BlobTools pipeline (version 1.1.1). BlobTools detects contaminants by searching with BLAST + (version 2.10.0) and evaluating the coverage of the sequence map against the assembled genome draft. Contigs identified as contaminants were removed using BlobTools, and contamination of the final genome sketch was assessed using CheckM pipeline with an independent method. This method uses marker genes specific to a given genomic lineage within a reference phylogenomic tree to estimate genome completeness and contamination. CheckM uses Prodigal (version 2.6.3) to identify genes in the genome draft and obtain the amino acid sequence of the predicted genes, HMMER (version 3.2.1) to identify the families of the identified proteins, and pplacer (version 1.1. alpha19) to determine the positions of these proteins in a pre-calculated phylogenetic tree. Taxonomic identification was performed using the Ribosomal Multilocus Sequence Typing online service (rMLST) (https://pubmlst.org/rmlst/) and searches against PubMLST using the mlst program (version 2.19.0) (https://github.com/tseemann/mlst).
Prokka annotation pipeline (version 1.14.6) was used for prediction and functional annotation of several gene categories in the genome by combining different tools: (1)  Absolute frequencies of the Enterococcus sp. were compared between fresh and dry samples (C). *indicates significant differences (Mann-Whitney test) between the groups, p<0.05. The data were tabulated in Microsoft Excel and examined using IBM SPSS Statistics for Windows, Version 21.0 (IBM Corp.) and GraphPad Prism version 7 (GraphPad Software, Inc., San Diego, CA, USA). The frequencies of the different findings were determined, and statistical differences were evaluated for CFUs and relative frequency of isolation. Normality and homogeneity of variance was calculated using the Shapiro-Wilk normality test and F test, respectively. After confirming the non-Gaussian distribution, a test for comparing independent groups was performed (Mann Whitney test).
Among the fecal samples, 29 showed Enterococcus growth at the following frequencies: E. dispar / E. durans / E. hirae (37.4%), E. faecium, E. gallinarum, E. casseliflavus (3.45%), and Enterococcus not identified at the species level (51.7%). [ Figure 1A]. The samples positive for E. faecium were from fresh stool. In addition, one sample from dry stool tested positive for S. aureus. A greater CFU count was observed in dry droppings than in fresh droppings (p<0.05) [ Figure 1B]. Species variation was not dependent on the type of sample (dry × fresh; p>0.05) [ Figure 1C].
S. aureus was susceptible to all tested antibiotics. In contrast, phenotypic resistance to gentamicin, erythromycin, and ciprofloxacin was observed in E. faecium, as well as susceptibility to vancomycin and linezolid.
The metrics for S. aureus and E. faecium genome assembly, including the total size, N50, and E-size (or auN), are shown in Table 1. Taxonomic identification by rMLST classified one of the isolates as S. aureus, with 55 identified alleles, all from S. aureus. MLST results indicated that the isolate belonged to S. aureus strain type 188 (ST188), with the following alleles: arcC (3), glpF (1), aroE (1), gmk (8), tpi (1), pta (1), and yqiL (1). Moreover, no alleles from other species were identified. No significant contamination was identified using either BlobTools or CheckM (Supplementary Table 1). Among the 1169 S. aureus markers present in the CheckM database, 1154 were found in a single copy, 10 were not found, and 5 were duplicates. This indicates good completeness and a low level of contamination.
CheckM evaluation of the genomes showed a good degree of completeness and a low rate of contamination (Supplementary Table 1). The Results of BUSCO quality analysis are presented in Supplementary Table 2. A summary of the annotations for Prokka is provided in Table 2. The annotations of antibiotic resistance genes are provided in Supplementary Table 3 for each strain.
Thus, medically important gram-positive cocci were isolated by C. livia dropping samples collected from the external environment of hospitals in this study. Compared to other staphylococci, S. aureus ST188 is highly adapted to different host types. It has shown an increasing prevalence of community-acquired infection on the Asian continent, infecting farm animals, adults, and pediatric patients [8][9][10] . However, we have not found any studies to date that relate S. aureus ST188 to the colonization of C. livia birds. In connection with healthcare, this strain has been previously described as causing bacteremia in a patient who underwent hospitalization 11 . Despite its phenotypic susceptibility to the tested antibiotics, the isolate demonstrated genetic determinants related to resistance according to molecular analyses. Discrepancies between genotypic and phenotypic results are frequently reported in the literature, and this phenotype/genotype dissociation may be a consequence of genetic alterations promoted by the environment or by other genetic factors; that is, the expression of a resistance phenotype by an agent may be dependent on a certain environmental context and/or on the existence of other genes or mutations 12 .
Enterococci are normally present in the intestinal microbiomes of humans and other animals and play a prominent role in hospital environment-related infections 13 . Enterococci present an important challenge for antimicrobial therapy owing to their intrinsic resistance to several drugs, as well as their ability to acquire resistance from other bacteria 14 .
Despite presenting susceptibility to vancomycin, E. faecium presented phenotypic resistance to gentamicin, erythromycin, and ciprofloxacin and was found to harbor resistance genes for variousantibiotics (Supplementary Table 3). The strain found in this study, E. faecium ST1054, was previously identified as a nosocomial agent causing urinary tract infection in a patient admitted to a neurosurgery ward in Portugal 15 .
The presence of domestic pigeons in the external areas of hospitals is a potential risk factor for the transmission of pathogens among patients, especially those undergoing intensive care. Therefore, it is evident that measures should be adopted to control pigeon populations near areas where health facilities are located.

ACKNOWLEDGMENTS
This work was supported by Fundação de Amparo à Pesquisa do Estado Minas Gerais (FAPEMIG), National Council for Scientific and Technological Development (CNPq) and Coordination for the Improvement of Higher Education Personnel (CAPES; finance